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Abstract 

We relate the structure functions of deep inelastic lepton-nucleon scatter- 
ing to current-current correlation functions in a Euclidean field theory de- 
pending on a parameter r. The r-dependent Hamiltonian of the theory is 
P° — (1 — r)P 3 , with P° the usual Hamiltonian and P 3 the third component 
of the momentum operator. We show that small XBj in the structure func- 
tions correspond to the small r limit of the effective theory. We argue that for 
r — ► there is a critical regime of the theory where simple scaling relations 
should hold. We show that in this framework Regge behaviour of the structure 
functions obtained with the hard pomeron ansatz corresponds to a scaling be- 
haviour of the matrix elements in the effective theory where the intercept of 
the hard pomeron appears as a critical index. Explicit expressions for various 
analytic continuations of structure functions and matrix elements are given 
as well as path integral representations for the matrix elements in the effec- 
tive theory. Our aim is to provide a framework for truely nonperturbative 
calculations of structure functions at small XBj for arbitrary Q 2 . 
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1 Introduction 



In this article we shall discuss the small XBj behaviour of the structure functions of 
deep inelastic lepton-nucleon scattering (DIS). The findings of the experiments HI 
and ZEUS at HERA (for recent summaries see [[[]]- [H) have brought this topic to the 
forefront of theoretical interest. Soon more data will come from HERA2. There are 
numerous suggestions for the theoretical description of the small Xsj behaviour of 
the structure functions, see for instance (|, || ||. Let us just mention a few of these 
approaches with representative references. 

As the first group of approaches let us mention the ones based on perturbative 
QCD, which allows one to derive evolution equations for the structure functions. 
It should be kept in mind that already in the derivation of these evolution equa- 
tions one has to make various assumptions and their practical use involves further 
approximations. 

Most popular and widely used is the DGLAP equation @ to calculate the evo- 
lution of the structure functions with Q 2 (see for instance 0, [Tl |3|). Improvements 
of the DGLAP method in fixed order in the strong coupling parameter a s , involving 
resummations to all orders in a s , have recently been proposed §|. Another time- 
honoured approach is based on the BFKL equation [|T0[ which is described in detail 
for instance in ||11|| . However, very large higher order corrections have been found 
in this approach fl2]| . Different recipes for dealing with this problem have been pro- 
posed |~3] . 

Other approaches make more assumptions and could be called QCD-based or 
-inspired models. Very popular at the moment are dipole models |14], Other 



approaches are based on the semiclassical approximation |RjL |6|, [T7|, [18j and on the 
colour glass condensate idea |L9|, [20]. 



Quite a different type of approach is based on Regge theory. It was shown in 
21, 22] that the small xbj behaviour of the structure functions can be well described 
using two pomerons, a hard and a soft one. Both pomerons are assumed to behave 
like simple Regge poles with linear trajectories. 

In this article we continue the investigations of the approach p3| where the be- 
haviour of the structure functions at small xbj is related to that of matrix elements 
in an effective Euclidean field theory. In this was explored for a model scalar 
field theory and it was argued that the limit Xsj — > corresponds to critical be- 
haviour in the effective theory. Here we extend these considerations to the case of 
QCD. The aim of this approach is to provide a framework where the small xbj be- 
haviour of the structure functions can be calculated from first principles using truely 
nonperturbative methods, for instance lattice methods. 

Our article is organised as follows. In section 2 we discuss kinematics, the reduced 
matrix elements free of kinematical singularities and the analytic continuation from 
the real to the imaginary z/-axis. Section 3 deals with the Deser-Gilbert-Sudarshan 
(DGS) representation which we use for further analytic continuations. In section 4 
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we discuss our effective Hamiltonians and Lagrangians for QCD, both in Minkowski 
and Euclidean space. Section 5 contains phenomenological applications and section 
6 our conclusions. 

2 Kinematics and analytic continuation in the v- 
plane 

In this section we recall first some definitions and results from DIS. The central ob- 
ject of our study is the forward virtual Compton scattering amplitude for nucleons. 
The absorptive part of this amplitude gives for real v the measurable cross sections 
of DIS. Our theoretical investigations will concentrate on the amplitude at imagi- 
nary values of u, from which nevertheless we can obtain information on these cross 
sections. 

2.1 Kinematics 

We study the forward virtual Compton scattering amplitude (Fig. 1) 

l*(q) + N(p)^ 1 *(q) + N(p), (2.1) 

where N stands for proton or neutron. We consider only spacelike virtual photons, 
q 2 = —Q 2 < and the amplitude averaged over the nucleon spin. The familiar 
Feynman amplitude for ( |2.1| ) is 

T^ipA) = ^u\ d ^ ^\ E (N(p)\T*J,(x)J u (0)\N(p)). (2.2) 

spins 

Here M is the nucleon mass, J /t (x) is the hadronic part of the electromagnetic 
current, and T* indicates the covariant version of the T product (see for instance 
||24j| ). It is understood that only the connected part of the matrix element is taken. 
All our conventions on kinematics follow 

In the following we shall, however, not work with ( |2.2|) but with the retarded 
amplitude 

T P(P^) = ^S d ' x ^ £ (N(p)\9(x°)[J,(x),J u (0)] cov \N(p)), (2.3) 

spins 

where we define 

0(a: o )[J„(x), J,(0)U = T*(J A4 (x)J I/ (0)) - J,(0)^(x). (2.4) 
The standard W^ u tensor of DIS is 

W^p, Q) = ^J d 4 xe^M, u (x,p), (2.5) 
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where 

M^(x,p) = \ £ (N(p)\J,(x)UO)\N(p)). (2.6) 

spins 

The expansion of Wnv in terms of the structure functions W\ 2 reads for pg > 0: 
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' 0"^ \Pu T- W 2\V, Q 



' M 2 \^ q 2 ; V q 2 , 

v = pq/M, Q 2 = -q 2 . (2.7) 

In the usual way we define the structure functions ^1,2,1, as 

F!(x Bj ,Q 2 ) = 2MW 1 {u,Q 2 ), 
F 2 (x Bj ,Q 2 ) = vW 2 (v,Q 2 ), 

AM 2 

F L (x Bj ,Q 2 ) = F 2 (x Bj , Q 2 ) - (1 + -Qz-x 2 Bj )~ 1 x B jFx(x B j, Q 2 ), 

and the transverse and longitudinal cross sections 

-^a T (v,Q 2 ) = W x {v,Q 2 ) 
Air z a 

' Fi(x Bj ,Q 2 ), 



2M 

K -a L (u : Q 2 ) = 1 ^±^W 2 (u,Q 2 )-W 1 (u,Q 2 



4vr 2 a " v ' Q 



1 , 4M 2 2 , „ , , z 
l + - F ^x%)F L (x Bj ,Q' 



2Mx B j y Q 1 
O 2 

K = v-£-. (2.9) 

2M v ; 



The expansions for T* v (|2.2j ) and T r ^ (|2.3| ) read 

Tjr (p, ff ) = (-^ + ^) if re v, g 2 ) 



1 



' M 2 J J T 2 ( 2 - 10 ) 

where now —00 < pq < 00. For Q 2 fixed, Tj{y,Q 2 ) and Tj ct (u,Q 2 ) (J = 1,2), are 
limits of functions Tj(i>,Q 2 ) analytic in the cut //-plane as shown in Fig. 2. The 
position of the nucleon poles is v = ±Q 2 /(2M), the cuts start at 

u = ±u t = ±[Q 2 + (M + m^) 2 - M 2 ]/(2M), (2.11) 
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where is the pion mass. We have for real v and j = 1,2 



Tf(u,Q 2 ) = limT>(l+ie),Q 2 ), (2.12) 



y2\ 



Tf{u,Q 2 ) = \im T^ + i^Q 2 ), (2.13) 
lmTf(u } Q 2 )=e(u)W j (u,Q 2 ) + e(-iy)W j (-u,Q 2 ), (2.14) 

imr; et (^g 2 ) = e^w^Q 2 ) -e{-v)w 3 {-v,Q 2 ). (2.15) 

For our purpose it is convenient to define in addition the following scalar ampli- 
tudes: 



v 2 + Q 2 



3T[>, Q 2 ) - _^LT 2 re >, Q 2 ), (2.16) 



T b re >,Q 2 ) = M~ 2 p^T^(P,q) 



fJ,U 

" 2 Q2 -Tr^Q 2 )+ f^^V T 2 rCt ^^')- (2.17) 



In a similar way we define for real v the functions T^ b and W a ^ and for complex 
v the functions T^iy, Q 2 ). The latter are analytic in the cut //-plane and have as 
boundary values for Im v — > +0 the functions T^^z/, Q 2 ). 

2.2 Analytic continuation in the z/-plane 

We are interested in the behaviour of Wi^i^yQ 2 ) for fixed Q 2 > and z/ — > oo. 
Instead of investigating the structure functions directly we shall study first the 
behaviour of the functions Ti i2 (^, Q 2 ) or equivalently T a ^{y, Q 2 ) for large imaginary 
u, that is for v = ir] with t] — > oo. Then we use the Phragmen-Lindelof theorem (see 
theorem 5.64 of |26j) to relate the behaviour of the amplitudes for large real and 



imaginary values of v. This has been discussed in detail in appendix A of [^3 



We follow now the same strategy as in [231]. We work in the rest system of the 



nucleon and choose the direction of q as third axis of the coordinate system: 



Using rotational symmetry and the same arguments as in p3f , we can represent the 
functions T aj h(z/, Q 2 ) for Im v > as follows: 
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1 /"CO poo i 

T a {y, Q 2 ) = - »,,■>, / dx° / dx 3 x 3 exp[ix°u - ix 3 Ju 2 + Q 2 



X 



I E<^(p)l(-^)^ )[J M (^es,^),^(0)]«»|^(p)>, (2.19) 



2 

spins 



j /"OO /"OO / 

Q 2 ) = - — . , =; / dx° / dzV exp[ix°u - ixV^ 2 + Q 2 
M\/v z + Q*Jo J-oo v 



x| E (N(p)\6(x°)[J (x 3 e 3 ,x ), J (0)U\N(p)}. (2.20) 

spins 

For v = %r\ and 77 > Q we get 

T a (zr/,g 2 ) = -^===^ dx° j_^dx 3 x 3 exp[-x° V + x 3 ^ V 2 -Q 2 } 

E(^)l(-^)^(^ )[^(^ 3 e3,x ),J,(0)] co ,|iV(p)), (2.21) 

spins 

T b {ir],Q 2 ) = My / J dx ° J_ oo dx 3 x 3 exp[-x°T] + x 3 ^r] 2 -Q 2 } 

x| E (N(p)\e(x°)[J (x 3 e 3 ,x ), J (0)}cov\N(p)). (2.22) 

spins 

The commutators in fl2.19|) -( p.22 ) vanish for \x 3 \ > x° and thus the effective inte- 
gration range is |x 3 | < x . This makes the integrands in (|2.19|) , ( |2.20|) exponentially 
damped for Im v > and represents the standard way of analytic continuaton of 
T^(u, Q 2 ) into the upper half i/-plane. However, the singularities on the lightcone 
|x 3 | = x° make it advisable to keep the x 3 -integrals in ( |2.19| )-( |2.22| ) as running from 
— 00 to +00. From T ab (irj,Q 2 ) we can derive T li2 (i77, Q 2 ) using ( |2.16| ), ( p. 17] ). 

In the next section we will derive another representation for T 12 (i^, Q 2 ) in which 
the integration path avoids the light-cone. 



3 The matrix element 

In this section we list general properties of the matrix element J\A VLV [x^p) ( [2.6|) and 
discuss its analyticity properties using the Deser-Gilbert-Sudarshan (DGS) repre- 
sentation WK, I2HI. 



3.1 General properties of M. llv {x, , p) 
Hermiticity of the electromagnetic current gives 

M„ w (x,p)* = M^(-x,p). (3.1) 
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Time reversal invariance of QCD requires 

M^(x,p)=M^(x',p'), (3.2) 



where x ^ = x^, p^ = p^. 
Current conservation implies 



d -M^(x,p) = 0, 



dx^' 

^M^(x,p) = 0. (3.3) 

For x 2 < the currents J M (x) and J u (0) in (|2.6| ) commute, which leads to 

M^ u (x,p) = M vfl (-x,p), for x 2 < 0. (3.4) 

The expansion of M.^ v {x^p) in terms of two causal scalar functions which are free 
from kinematical singularities and where ( |3.1|) to ( |3.4|) are satisfied automatically 



[29| reads 

M^x^p) = [g^D - d^Miixp, x 2 ) 

+ [PvPv u ~ (pd)(ppd u + Pvdp) + g^(pd) 2 ]M 2 (xp,x 2 ). (3.5) 

The invariant functions Aij (j = 1, 2) satisfy: 

Mj(xp, x 2 ) = M*(-xp, x 2 ) (3.6) 

and 

Mj(xp,x 2 ) = Mj{-xp,x 2 ) for x 2 < 0. (3.7) 

3.2 Representations of amplitudes in terms of 

Inserting ( |3.5|) into (|2.5|) we obtain the structure functions Wi j2 of (|2.7|) in terms of 
the Fourier transforms of A4i$'- 

W^Q 2 ) = J d'xe^M^x 2 ) + d 4 xe^M 2 (x Pl x 2 ), (3.8) 

W 2 (u,Q 2 ) = ( d A xe iqx M 2 {x P) x 2 ). (3.9) 

Note that in (|3.9|) the kinematic zero of W 2 {v, Q 2 ) at Q 2 = is made explicit, as it 
should be. For the cross sections (|2.9| ) we get: 

-^-a T (v,Q 2 ) = — L- fd 4 xe iqx [M 2 u 2 M 2 (xp,x 2 )-Q 2 M 1 (xp,x 2 )},(3.10) 
K , ^ Q 



4tt 2 



a 



cr L (u,Q 2 ) = ^j^xe^[M 2 M 2 {T P ,^)+M x {x Pi x l )\. (3.11) 
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Inserting ( |3.5| ) into Q2.3| ) and using ( 2.10|) and rotational symmetry in the nucleon 
rest frame, we get 



1 



poo poo _ i 

/ dx° / dx 3 x 3 exp[ix°is - ix 3 J v 2 + Q< 

JO J-oo v 



M^/v 2 + Q 2 Jo 
x{-Q 2 [Mi(xp,x 2 ) - Ml(xp,x 2 )] 

+M 2 v 2 [M 2 (xp, x 2 ) - M* 2 (xp, x 2 )]}, (3.12) 
MQ 2 



roc roo _ y 

/ dx° / dx 3 x 3 exp[ix°z/ — ix 3 J v 2 + Q 2 ] 

JO J-oo v 



where 



xp = x u M , x 2 = (x ) 2 - (x 3 ) 2 



(3.13) 



(3.14) 



In the following we will make a change of variables as in [23| and replace (x°, x 3 ) 
by (x°, r), where 

(3.15) 



= 1 - — 

r x ' 



We define 



Mj{x°, r) : = Mj(x°M, (x°) 2 r(2 - r)) 
(J = l,2). 



With this we get 

Tl et (u,Q 2 ) 



1 r°° 
- MVu 2 + Q 2 L dr{1 - r) 

x y°° rfx°(x ) 2 exp[ix°(z/ - (1 - r)^v 2 + Q 2 

x{-Q a [^x°,r)-(^rOr ,r))*] 
+ MV[A^ 2 -(x°,r) - (>l 2 -(x ,r))*]}, 



(3.16) 



(3.17) 



T 2 r >,Q 2 ) 



MQ 5 



rfr(l — r) 



y/u 2 + q 2 j-oo 

/•CO / 

x jf dx (x°) 2 exp[ix°(z/-(l-r) v /z/ 2 + g 2 
x[M^(x ,r)-(M^x°,r)n 



(3.18) 



3.3 DGS representations for 

For the causal functions M.\^{xp, x 2 ) defined in (|3.5| ) we can write down DGS rep- 



resentations |27|, g8| . The behaviour of the structure functions for small XBj as mea- 
sured in experiments (see chapter 5 below) suggests that the DGS representation 
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for A4 2 needs no subtraction, whereas .Mi needs subtractions. Thus we write 

M 2 {xp 1 x 2 ) = ds d(f 2 (s,()exp(i(px)-A + (x,s), (3.19) 
Jo J-i i 



where 



(2tt) 3 
1 



(fk e- zkx 6(k u )5(k 

1/2 



An 2 \—x 2 + iex° 



K x Js(-x 2 + iex°) 



(3.20) 



is one of the usual invariant functions given in terms of the modified Bessel function 
K\. In ( |3.19|) f 2 is a function for which ( |3.6|) and ( |3.7|) require 



/2(a,C) = /a(a,-C)=^(a,0- 
Thus, ( |3 . 1 9| ) can also be written as 

/■oo />1 I 

Jvl 2 (xp,x 2 ) = / / dC,f 2 {s, C) cos((jxc)-A + (:r, s) 
Jo J-l z 

For M.\ we write down a subtracted DGS representation: 



(3.2i; 



(3.22) 



J\4i(xp, x 2 ) 
where 



ds\fi°\s)+ I da\(s,Oicos((px) - 1} \ -A + (x,s) (3.23) 



A w « = (/{%))*, 

/ 1 (a,0 = /i(a,-0 = /r(« > 0- 



(3.24) 



In P7| , P8| the DGS representation is primarily discussed for the matrix element 
of the commutator and the T-product of two currents. From the DGS representation 
of the T-product we get immediately the representation for the ordinary products 
of currents used here. 

Let us now switch to the matrix elements A4~(x°,r) fl3.16| ) to discuss their 
behaviour in r for fixed x°. For M. 2 we get from ( |3.19| )-( p.22| ) for a; > 0: 



M 2 (x°,r) = M 2 (x°,r,e) 

poo rl 

= ds d(f 2 (s,()cos((Mx°) 
Jo J-l 



1 s 1 ' 2 



X 



An 2 x' 



_ (_ r ( 2 - r ) + ley 1 ' 2 Kt (s 1/2 x Q {-r{2 - r) + ief 



/2 



(3.25) 
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For clarity we keep in the following discussion the e-parameter (e > 0) as an explicit 
argument. A similar expression is obtained for hA\ from ( |3.23|) . 



We will now make the assumption that the weight functions fj(s,Q are suffi- 
ciently well behaved such that the representation ( |3.25| ) for Ai2{x°,r,e) and the 



analogous one for .Mi allow an analytic continuation in r at fixed x° and e. The 
singularity structure of A4j(x°, r, e) in the r-plane can then be read off from ( |3.25| ). 
There are two branch points, at r = 1 =p y/1 — ie, giving for e«l, 

r ~ — e , and r « 2 — -e. (3.26) 

We draw the associated cuts to the right (see Fig. 3a). The cut structure of 

(M 3 (x°,r*,e)Y =Mj{x°,r,-e) (3.27) 



is shown in Fig. 3b. Looking now at (|3.17 ) and ( |3.18| ) we see that the r-integrations 
run on the real r-axis from r = — oo to r = +oo over Aij(x°, r, e) and from r = oo 
back to r = — oo over A4j(x°,r, — e). Deforming the contours slightly we can obtain 
T[% as integrals over the analytic functions 

Mj(x°, r) := Mj(x°, r, 0) (3.28) 

along a curve C (see Fig. 4). The curve C comes from the second sheet, moves onto 
the first sheet between r = and 2, circles r = and moves again onto another 
sheet between r = and 2. We get 

t; ^- q2) = -M^rwL ir[l - r) 

x dx°{x°) 2 exp[ix°{u - (1 - r)^Ju 2 + Q 2 )] 

x {-g 2 A^i(s°, r) + M 2 v 2 M< 1 {x\ r)}, (3.29) 

x dx°{x ) 2 exp[ix°{u - (1 - r)^u 2 + Q 2 )] M 2 {x°,r). 

(3.30) 



o 



Note that the matrix elements M.^ (x°,r) of p. 16] ) for < r < 2 are the limits of 
the analytic functions A4j(x°,r) below the cut: 

Mj(x°,r) = lim Mj(x°,r-ie), (0 < r < 2). (3.31) 

J £— »+0 

We can now use ( 3.29 ), ( 3.30 ) to perform the analytic continuation in v to the 



upper half z/-plane in an alternative way to that of section 2.2. We move the curve 
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C arbitrarily close to the positive real r-axis, r > 0, and use the inequalities (A.l) of 
p3fl to show that the exponentials in ( |3.29| ) and ( |3.30|) are damped for Im v > 0. We 



are interested in Ti^ for large imaginary v. With v = irj, i] > Q, the exponentials 
can be written as 



exp[ix (u - (1 - r)\Jv 2 + Q 2 )\\ v=ir) = exp[-x°/x(r], r)], 



x(rj, r) — [r) — (1 — r)\Jrf — Q 2 ] 



2i-i 



1 1 + f 

r] r + f ' 



277 2 i _ ^ + /i 



(3.32) 
(3.33) 

(3.34) 



With (|3.32|) we get for the analytically continued functions 7\ 2 from ( |3.29| ) and 

dH3) 



Ti(^,g 2 ) 



x exp[-a; /^(^, r)] [g 2 A^i(s°, r) + M 2 rj 2 M 2 (x°, r)], (3.35) 



T 2 (^,Q 2 



x exp[— x°/x(r), r)] A4.2{x [i ^r). 



(3.36) 



We note that for given Q 2 and r\ > Q the characteristic damping length x(j],r) 
is positive for 

Re r > -f. (3.37) 

We will thus take in the following the curve C in ( |3.35| ) and ( |3.36| ) to be to the right 
of — f as shown in Fig. 4. In this way we get integral representations of T 12 (irj, Q 2 ) 
where the integration path avoids the singularities of the integrand at r = 0, 2, which 
correspond to the light-cone. 

The representations ( |3.35| ) and ( |3.36| ) will be used as basis for the discussion in 
the following sections. 

Let us now estimate the relevant integration range in r in ( 3.35Q and ( 3.36Q for 
f] — > oo. We see from ( |3.33| ) that for fixed r we have 

_1 for r] 



x(r], r) 



V 



oo. 



(3.38) 



Then the factors exp(— x°/x(rj, r)) will suppress such contributions to the integrals 
( |3.35| ) and (|3.36|) . Now we can keep the curve C in Fig. 4 at finite fixed values of r 
except for the region between r = and r = — f where C has to cross the negative 
real axis. Taking typical value r = — r/2 we get 



x(r), -\r) 



2(1 + f) 



2rj 

Q 2 f ° r V 



oo. 



(3.39) 
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Thus, for i] — > oo this region, where x(r],r) becomes very large, will give the main 
contribution to the integrals ( |3.35|) and ( |3.36| ). The behaviour of Ti^iirj.Q 2 ) for 



77 — ► oo is therefore expected to be governed by the behaviour of Aii t 2(x°,r) for 
small |r| and large x°. 

4 Effective Hamiltonians and Lagrangians 

In this section we shall express the matrix elements M.^ (|3.5| ) as functional inte- 
grals with effective Lagrangians containing r ( |3.15| ) as parameter. The procedure is 
analogous to the one of ]23| but the vector nature of the electromagnetic current 
causes some complications. We work in this section with the Lagrangian of QCD in 
a general covariant gauge. The case of the temporal gauge is treated in appendix C. 

4.1 Matrix elements in Minkowski and Euclidean space 

We start here with the following matrix elements in the nucleon rest frame, supposing 
always x 3 = x°(l — r) , x° > 0: 

M-{x\r) = \ £ (N(p)\(-gnMx 3 e 3 ,x°)J u (0)\N(p)}, (4.1) 

spins 

M b (x°,r) = \ J2 (N(p)\J (x 3 e 3 ,x°)J (0)\N(p)). (4.2) 



spins 



The relation of A4 ab (x°,r) and Ai 12 {x°,r) ( |3.16|) is given in appendix A. From 
translational invariance we have 

J u (x 3 e 3 , x°) = exp(ix°H r )J^(0) exp(-ix°H r ), (4.3) 

where 

H r = P°-{l- r)P 3 (4.4) 

and P°, P 3 are the energy and third component of the momentum operator. Note 
that H r is positive semidefinite for < r < 2. From ( |4.3| ) we get 

M-(x°,r) = lY.( N (p)\(-9neM™°Hr)Ju(0)exp(-ix°H r )J u (0)\N(p)), 

spins 

(4.5) 



M^(x°,r) = | ^(iV(p)|exp^ J ff r )J (0)exp(-^ i7 r )J (0)|7V(p)). 



spins 



(4.6) 



We see that in the effective theory Ai ab (x°,r) are correlation functions of two cur- 
rents at purely timelike separation x° . 
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We will now keep r fixed with < r < 2 where H r is positive semidefinite. This 
allows us to continue A4~ b (x°, r) analytically into the lower half x°-plane and in 
particular to x° = — iX 4 with X 4 > 0: 

M-(-iX,,r) = \ Y, <JV(p)l(-0 exp (X A H r ) J„(0) exp (-X 4 H r ) J v (0)\N(p)), 

spins 

(4.7) 

M^(-iX 4 ,r) = \ ^(iV(p)|exp(X 4 if r )Jo(0)exp(-X 4 ^ r )Jo(0)|iV(p)). 

spins 

(4.8) 



The matrix elements (|4.7|) , ( |4.8|) are correlation functions of two currents at purely 
timelike distance X 4 in a Euclidean field theory with H r as Euclidean Hamiltonian. 

In the following we shall write the matrix elements (|4.5| ) and (|4.6| ) as path in- 
tegrals in an effective, r-dependent, Minkoswkian theory and ( |4.7p and ( [4.8| ) in an 
effective Euclidean theory. 



4.2 Effective Lagrangian and path integral in Minkowski 
space 

We start with the Lagrangian of QCD in a general covariant gauge 

C = -\G%G^-\r\d^Gl){d v Gl)+Y,^D,-m q )q 

+ (d^ a )(d a 4> a ) ~ gfaUd^ a )G b u <P c . (4.9) 

Here q denotes the quark fields, 4> a the Fadeev- Popov fields, G° the gluon potentials 
and G" the gluon field strength tensor, 

G% = d,G a v - d u Gl - gf ahc G\Gl. (4.10) 

The coupling constant is g, the quark masses are m q , and is the covariant deriva- 
tive 

D a q={d a + igGl\\ a )q. (4.11) 
All quantities in (|4.9| ) are the unrenormalised ones. 



It is now a straightforward exercise to derive from ( f4.9|) the canonical momenta II, 
the Hamiltonian density, the Hamiltonian P°, the third component of the momentum 
P 3 and the effective Hamiltonian H r ( [4.4| ). The details are given in appendix B. Here 
we list the result: 

H r = I d 3 xH r , (4.12) 
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7~C r = — ^£II(3aoII(3aO 

-n Ga0 {djG^ - (1 - r)d 3 G a0 ) 

+ ^H G ajU.Qaj 

+U G a i {-d j G a0 + gf abc G b0 G^ + (1 - r)d 3 G a n 
+\G%G a > k 

+ E 7 J % +1 7°(1 - r) H 3 +^7 M G^A a + m q )q 
q 

+n^n^ + n^(^/ o6c G 6 V c + (i - r)<9 3 a ) 

+ (l-r)(9 3 a )n^ + (^ a )^ a 

+gf abc {d j r)G b ^ c . (4.13) 

Here and in the following Latin indices j, k, ... run from 1 to 3. From the effective 
Hamiltonian density ( |4.13| ) we get in the standard way the effective Lagrangian 
density: 

C r = -±Gl v G a » v - (1 - r)G a ° 3 d 3 G aj + 1(1 - r) 2 (d 3 G aj )d 3 G aj 
-\C\d^G a » - (1 - r)d 3 G a0 )(d u G av - (1 - r)d 3 G a0 ) 

+ £ ^ B M -m, - (1 - r)f 7 ° d 3 )<7 

+(«9> a )<9 M a - gf abc (d^ a )G b ^ c 
-(1-r) ((9 3 a )0 ffl + ^ 3 a ) 
+ (l-r) (<9 3 a ) ^ afec G w c 

+ (l-r) 2 to)9 3 r- (4-14) 



Here and in the following we use the term "effective" for the r-dependent theory. 
We emphasize that our procedure does not imply approximations like integrating 
out some modes etc. The r-dependent theory is as good as the original one for the 
calculation of our matrix elements and we hope that it will be more effective for the 
study of the small XBj limit of the structure functions. 

We can now give the path integral representation of Ai~ b (x°, r) (|4.5| ), fl4.6|) in the 
effective theory described by ( |4.14| ). The procedure is analogous to the one described 
in detail in f23|. Let iPn(x) be an interpolating field operator for the nucleon with 
normalisation such that 

(0\^ N (x)\N(p,s)) = e^u.ip), (s = ±1/2). (4.15) 

According to the LSZ formalism [3(| we define 

A,(p,t)= I d 3 xu s (p)e ip *j°ilj N (x). (4.16) 

Jx°=t 
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In the sense of the weak limit we have 



lim A s (p,t) = A?(p), 

t — >—oo 

lim A s (p,t) = A^(p), (4.17) 

t — >+oo 

where A l f-{p) and A° ut (p) are the annihilation operators for incoming and outgoing 
nucleons. In the standard way we obtain now for the matrix elements 

M~(x°, r) = t lim^ exp[iM(t f - U)} 

ty — > + oo 



X 2 



:Z 1 J V(G, q, q, 0, 0) exp[i J d 4 xC r (x)] 

\ E as(p,t f )(-gnMO, x )j u (0)al(p, U), (4.18) 

s 

Z = Jv(G,q,q,<p,<i))exp{i J d A x£ r (x)}, (4.19) 
Mb(x°,r) = lim exp[iM(tf-ti)] 



t j- — * + oo 

X 



Z 1 J V(G, q, q, (ft, (j)) exp[i J d A xC r (x)] 
x \ E <P> tfMO, x°)j (0)al(p, U). (4.20) 



Here p is always the nucleon momentum in the rest system (see ( [2.18] )) and a s , o^ s ,j^ 
are obtained by replacing in A S ,A\, J M the quark field operators by the corresponding 
Grassmann variables. Of course, also in C r (x) quark, ghost and gluon field operators 
have to be replaced by Grassmann variables and classical gluon fields, respectively. 
With ( |4.18 )-( 4.20|) we have the path integral representation of the matrix elements 



Ai ab as timelike correlation functions of the currents in the effective, r-dependent 
Minkowskian theory. 

4.3 Effective Lagrangian and path integral in Euclidean 
space 



Here we derive the path integral representation for the matrix elements 
in the Euclidean effective theory. Points in Euclidean space are denoted by X = 
(X, X4). Our Euclidean 7-matrices are 

lEj = ~il\ (j = 1,2,3), 

7E4 = 7°- (4.21) 

We perform now a rotation to Euclidean space starting from the effective Hamilto- 
nian (|4.12p , ( 4.13|) . The details are given in appendix D. The result for the effective 
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Lagrangian density in Euclidean space is: 

C E , r = IrH^ + *(1 - r)d 3 G a E4 ][duG a Eu + z(l - r)d 3 G^ 4 ] 
+±[G| 4j + i(l - r)«9 3 G^][G^ + i(l - r)d 3 G a Ej ] 

I 1 /^a /^ya 

+ 53 ( 57Ka» <9 m +57^4(1 - r) <9 3 +ig-iE^G a Eil \\ a + m q ) q E 



+(dS a E + i(l - r)d^ E ){d A r E ~ 9fabcG b m( p E + i(l - r)«9 3 <^), 
+(fi^|)a^J - gUdd^G^. (4.22) 

For the following it is convenient to split £#,r into the quadratic and the inter- 
action term: 

£e,r = £>E,r + ^E^r) (4.23) 

+|[a 4 G^ - djG% A + z(l - r^G^RG^. - + i(l - r)<9 3 G^.] 

+ 2^ ^[2^ 9/i +57134(1 - r ) ^3 +mg]g E 
+(d^ a E + z(l - r)d 3 4> E ){d^ a E + i(l - O^) 

+(W(^), (4.24) 
—\gfabc{djG a Ek — dkG Ej )G b Ej G c Ek 

1 1 2 /" /" /""fC y-fil 7 /^YC' 

"T4i/ J abc} ab' c'^> Ev^ Ep^ Ev 

+ig J2 QElE^G Epi \\ a qE 
1 

-gfa bc [(d,r E )G E ^ + z(l - r)(«9 3 <^)G^ 4 ]0 c B . (4.25) 



With the same procedure as in j£J we get the following path integral representation 
for the matrix elements -M~ fe ( — iX&, r) of ( |4.7|) and (|4.8| ): 

.M~(-Or 4 ,r) = |^ lim exp[(r f - Ti)M] Z E X V(G E ,qE,qE,4>E,4>E) 

Ty? — > + 00 

xa S)B (p,r/)(-l)jB At (0,X4)j S/1 (0)4(p,ri)exp(-5B, r ), (4.26) 
.M^(-zXi,r) = |E lim ex P[( r / - T i)M] Z E l V(G E ,q E ,q E ,(f) E ,(j) E ) 

O 7- j. — too 

xa s , E (p, Tf)jE,4(0, X 4 )j E4 (0)aJ(p, n) exp(-S , jB)r ). (4.27) 
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Here 



Z E = J D(G E ,qE,qE,(f>E,<f>E) exp(-S E>r ), 
Se,t — J d A X C Er (X), 



(4.28) 
(4.29) 



G E ,...,(p E are the integration variables, j EfJ , represent the components of the elec- 
tromagnetic current and a S) a\ the nucleon state; see appendix D. 

From (|4.23[ )- fl4~2"5| ) we see that the Euclidean, r-dependent action S E , r has, in gen- 
eral, imaginary parts. To discuss the formal convergence properties of the integrals 
over the gluon-potential variables we split S E>r into the quadratic and interaction 
parts 

S E , r = 4°l + ^ (4.30) 



q(0) 
°E,r 



d 4 XC {0) 



E,ri 



E,r- 



(4.31) 
(4.32) 



After some partial integrations we get 
S E j r - J d 4 X^G a E ^5 ab 



- 5^d x d x + (1 - C l )d,d u - 2i(l - r)6^d 3 d 4 

- r)(l - £- l ){5^d v dz + 8 u4 d„d 3 ) 
+ (1 - r)H, v d 3 d 3 - (1 - r) 2 (l - C l )5^5 ui d 3 d 3 ]G h vE 

+ QEbEndp + i(l - r)-i Ei d 3 + m q ]q E 
q 

+r E 5 ab (-d x d x - 2i(l - r)d 3 d 4 + (1 - r) 2 d 3 d 3 )<p b E } . (4.33) 

It is particularly convenient to use the Feynman-'t Hooft gauge, that is to set £ = 1 
in the following. With this we find 

4°lk=i = / d 4 x{|GVU<W-dAdA - 2i(l - r)«9 3 <9 4 + (1 - r) 2 9 3 9 3 ]^ 

+r E S ab [-d x d x - 2i(l - r)<9 3 <9 4 + (1 - r) 2 a 3 a 3 ]4} . 
With the Fourier transformations 

(Ge,(K)T 
q E (X) 

<Pe(X) 



r Eu 



(4.34) 



^ e G E(1 {K), 



(27T) 



— vre qE[tt), 



d 4 K JKXJa 



(2^) 4 



d 4 K -iKX^g 
(2tt) 4 C 
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a (iT) 



(4.35) 



we get 

[K 2 + 2i(l - r)K 3 K 4 - (1 - r) 2 K 2 3 ]G b Eu (K) 

+ Yl ^E^VlE^K^ - (1 - r)>y E4 K 3 + m q ]q E (K) 

+J a (K)5 ab [K 2 + 2i(l - r)K 3 K 4 - (1 - r) 2 K 2 ]4> b (K)} . (4.36) 

Note that 

K 2 + 2i(l - r)K 3 K A - (1 - r) 2 ^ 2 

= K 2 + K\ + r(2 - r)K\ + iT 2 + 2i(l - r)K 3 K 4 (4.37) 

has a positive real part for < r < 2 and K ^ 0. 

We discuss now the formal convergence properties of the gluonic integrations in 
( |4.26| )- (f08| ), setting £ = 1. Consider a given gluon potential G Efl (X) together with 
\G E (X) for all A > 0, to investigate how the integrand in the functional integral 
behaves for large potentials. Suppose first that 

U c G b E ^X)G Ev {X) $ 0. (4.38) 

Then the term in C Er ( |4.25| ) quartic in the gluon potentials ensures 

exp[-S E A XG )} oc e~ A4ci for A -> oo (4.39) 

with Ci > 0. On the other hand, if 

f abc G b Ell {X)G Eu {X) = (4.40) 

for all X, the positivity of the real part of the quadratic term of the action ( [4.36D 
ensures that 

exp[-Sj5 )r (AG)| €= i] oc e- x2c2 for A ^ oo, (4.41) 

with C2 > 0. Thus in any case the integrand of the functional integrals ([4.26|) -( f4.28| ) 
is damped exponentially for large gluon potentials. This should make the integrals 
well behaved after introducing some regularisation procedure, for instance a lattice 
regularisation. 

4.4 Propagators in Euclidean space and perturbation ex- 
pansion 

In this section we discuss first the lowest-order propagators in Euclidean space for 
the gauge choice £ = 1. The basic Green's function A E (X, m 2 , r) for mass m in the 
r-dependent theory is defined through 

[ — d\d\ - 2i(l - r)d 3 d A + (1 - rfd 3 d 3 + m 2 } A E (X, m 2 , r) = 5\X). (4.42) 
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The solution of ( |4.42|) as given in [53 



is 



A E (X,m 2 ,r) = J^-^ KX (K T A r K + 



m 



2\~1 



in 



471-2 



d 4 K ... 

(x T A; 1 x)- 1 / 2 K 1 (m(x T A; 1 x) 1 / 2 ) 



(4.43) 



where K\ is the modified Bessel function of order 1 and 



f 1 \ 

1 

r(2 - r) i(l - r) 

\ z(l - r) 1 / 



(4.44) 



At 1 



For m 2 = we get 
A £ (X,0,r) 



/ 1 

10 

1 -i(l-r) 

V -r) r(2-r) / 



(X^T 1 ^)- 1 



(4.45) 



471-2 



^(X 2 + Xf + Xj + r(2 - r)X 2 - 2i(l - r)X 3 X 4 )" 1 . 



(4.46) 

From ( 4.34 ) we get now easily the expressions for the propagators in lowest order. 
The gluon propagator is 



A i ^(X,r) = 5 ab 5 IIU A E (X,0,r), 
the quark propagator is 

S E \x,m q ,r) = [-Je^ - i{\ -r)^ E4 d^ + m q ]A E (X, m 2 , r), 
and the ghost propagator is 

AP ab (X,r) = 6 ab A E (X,0,r). 



(4.47) 



(4.48) 



(4.49) 



All these propagators ( |4.43| ) , ( |4.47| )-( |4T4^| ) have the property that for r — > they 
fall off more and more slowly in X4 direction, whereas their fall-off in X±,X2 and 
X 3 directions is independent of r. For the case m 7^ this is discussed in section 
5 of [53J. For the propagator of a quark of mass m q the correlation length in the 
directions Xj(j = 1, 2, 3) is 

4 0) = m~\ (4.50) 
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and in the direction X4 



Zf = [mjr(2-r)]- 1 / 2 



4 0) (2r)" 1/2 forr-^0. (4.51) 



For the massless case ((4.46|) there is of course no genuine correlation length, but for 
r — > the fall-off in X 4 direction is again slower by a factor (2r) -1 / 2 compared to 
the Xl,2,3 directions. 

We can now use the propagators ( }4.47|) -( f4.49|) to construct the unrenormalised 
perturbation series for the Green's functions in the theory described by Ce,t Q4.22| ). 



Since the Fourier transforms of the propagators (|4.43|) , (|4.47|) - (|4.49| ) have no un- 



wanted poles in momentum space, the convergence properties of Feynman integrals 
should be the same for arbitrary r as for the standard case r = 1. Thus we conclude 
that also the construction of the renormalised perturbation series should work for 
the CE,r theory as for the standard case. Of course, the purpose of our paper is 
not to advocate perturbative calculations starting with the Lagrangian density Ce,v 
(|4.23|) but to provide a framework for nonperturbative calculations. 



5 Phenomenological applications 

In this section we study the naive parton model from our point of view and make 
some speculations concerning a possible critical behaviour of the full theory with 
interaction for r — > 0. 



5.1 The naive parton model 

The central assumption of the naive parton model is that of free field or canonical 



behaviour of the product of currents near the light cone. In detail one assumes |29 
for Mi 2 of 



4tt 2 

Mi(xp, x 2 ) ~ hi(xp) A + (x, m 2 ) 

1 

~ hi(xp)(— x 2 + ie(xp))^ 1 , (5.1) 

M 2 (xp,x 2 ) ~ -h 2 (xp)8n 2 ds5'(s - m 2 )-A + (x,s) 

Jo 1 

~ \h 2 (xp) ln[m 2 (— x 2 + ie(xp))}, (5.2) 

where ~ indicates that only the leading term for x 2 — > is considered and A + is 

defined in ( [3.20|) . The parameter m represents a hadronic mass scale. Its precise 
value does not matter except that we will assume 

< m < M. (5.3) 
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According to the DGS representations ( |3.22| ), ( |3.23| ) the functions hi )2 (xp) can be 
represented as 

h(xp) = f 1 d(icos((px) - 1]^(C) + hf\ (5.4) 



where 



h 2 (xp) = d(cos((px)h 2 ((), (5.5) 



k(0 = k(-C) = h*(C), (7 = 1,2), 
h (o) = h f)* = CQnst ( 56 ) 

It is now an easy exercise to calculate the structure functions W\ t2 inserting the 
parton model ansatz (|5.1|) - (|5.6|) into ( |3.8|) , (|3.9| ). The result, expressed in terms of 
F 2 r of (pi), is 



F 2 {x Bj ,Q 2 ) = z4U 2 (z/,Q 2 ) ~27r 2 x Bj ti 2 (x Bj ), (5.7) 

2 ; 



F L (x B j,Q 2 ) ~ 4ir 2 x Bj hi(x B j), (5. 



where now ~ indicates the leading term for Q 2 — > oo. Of course, we get Bjorken 
scaling. From the well-known relations of the parton model we get the physical 
interpretation of the functions hi and h 2 as 

^%(0 = E^(o + E^.(o, 

3 j 

< C < 1, (5.9) 

where ej,Nj(() are the charges and distribution functions for the spin 1/2 partons 
and ij, Nj{() for the spin partons. 

Using an approximate parton model description, we conclude from experiment 
0-|| and the fits in [^] that the growth of F 2 for x B j — > is at most as 

F 2 (x Bj ,Q 2 ) ex (x Bj y a (5.10) 

with < a < 1. This implies from Q5.7p that at most 

K(o « icr~\ 

h 2 (Q oc |C|- a , (5.11) 

for ( — > 0. The structure function F^ is not well measured at high Q 2 . All indications 
are that Fl is small compared to F 2 . Thus a conservative estimate is that also Fl 
grows at most as 

F L (x Bj ,Q 2 )(x(x Bj )- a (5.12) 
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for xbj — > 0. From 



this implies that at most 

ho oc icr~ 2 



(5.13) 



for ( — > 0. We see that with ( |5.11| ) and ( |5.13| ) the integrals ( |5.5| ) and ( |5.4j ) are 
perfectly convergent for ( = 0. On the other hand, an unsubtracted integral as for 
h 2 would not be convergent for h\. These findings are our motivation to write the 
DGS representations (|3.22|) without, and ( |3.23j ) with subtractions. 

The purpose of our discussion of the parton model is to see how one obtains in 
this case a relation between the behaviour of the structure functions at large z/, that 
is small x B j, of the amplitudes Ti 2 (|3.35| ), (|3.36| ) at large r], that is large imaginary 
z/, and of the matrix elements M.\ )2 (|3.28|) and M.~ h ( [4.1|) , ( |4.2|) at small r. For 
illustration we consider only W 2 ,T 2 and j\A 2 and make the following simple ansatz 



with 



0<a<l, b > 1 , A>0. 



(5.14) 
(5.15) 



We get then for < ( < 1 
2vr 2 ^(C) 



A[i-c+ k a ai-c a )]c a - 1 (i-o b - 1 

AC 11 - 1 for small ( 



and from Q5.71) for small XBj 

F 2 (x Bj ,Q 2 ) ~ A(x B j)' 
which implies for large v at fixed Q 2 : 



W 2 (u,Q 2 



<?t{v, Q 2 ) + Q 2 



A 



(2MuY 

\Q 2 J ' 

, 2 A (IMvX 

4n v M 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



On the other hand, we have argued at the end of section 3 that the behaviour of 
T 2 (i?7,Q 2 ) for 77 — > 00 should be governed by the behaviour of A^^ , r) ( |3.28| ) for 
small \r\ and large x°. With the ansatz (|5.2|) we get 



A^ 2 (a; ,r) ~ \h 2 (x Q M) ln(-r) 

for small |r|, where from ( [5.5|) and ( ^.14j ) : 

A 



h 2 (x°M) 



7r 2 a Jo 



c2Ccos(0r M)(C- a -l)(l-C) b . 



(5.20) 



(5.21) 



22 



To see the large x° behaviour of h 2 (x°M), we write it as follows: 

A r°° 

h 2 (x°M) = — / dCcos(Ca; M)C" a e- c 

7T 2 a Jo 

/ dCcos((x°M)i(C a - 1)(1 - 0^(1 - - C a e- C ). 
n z a Jo 

(5.22) 

The second integral on the r.h.s. of ( |5.22|) vanishes faster than l/x° for x° — > oo, as 
we find from a simple application of the Riemann-Lebesgue lemma. The first integral 
then gives the leading behaviour for large x°: 

h 2 (x°M) ^-r(l - a) sm(la)(x M) a ~\ (5.23) 

Inserting Q5.20| ), ( |5.23| ) into (|3.36|) gives for large rj: 

2 iMQ 2 AM a ~ l . 7T 

T 2 (ir),Q ) ~ r(l - a) sm(-a) 



f roc 

x dr ln(— r) / dx°(x°) 1+a exp(— x°/x(r], r)) 
Jc Jo 



iMQ 2 A 



M a ~T(2 + a)r(l - a) sin(fa) 



x dr ln(-r)(r + f)~ 2 ~ a 
Jc 

■ (5 - 24) 



(5.25) 



With analytic continuation we find for large \v\ 

T 2 (u,Q 2 ) ~ A(cos(fa))- 1 [sin(fa) +^cos(fa)]i 

and for large real v 

W 2 (v, Q 2 ) = Im T 2 {v + ze, Q 2 ) ~ A 1 - . (5.26) 

Of course, this is in perfect agreement with ( |5.18| ). 

We have thus seen in this simple example that the power behaviour ( |5.17| ) of F 2 
corresponds to a behaviour 

M 2 (x°, r) oc (x°M) a ^ 1 ln(-r) (5.27) 

for small \r\ and large x° (see Q5.20D , ( |5.23|) ). If we assume that A4i(x°,r) has a 
similar behaviour, as is for instance true for vanishing <jl (see ( |3.11| )), we get from 
(|5.27| ) and ( |A.1|) , (|A.2|) for the current correlation functions ( |4.7|) , Q4.8|) : 

M- b (-iX 4 ,r) oc (-iX 4 M) a - 3 (-r)- 2 (5.28) 

for small \r\ and large X 4 . 



23 



5.2 Regge behaviour 

We have seen in section 4.4 that the free propagators in the r-dependent Euclidean 
theory develop in X 4 direction a large correlation length oc \jyjr for r — > 0. Let 
us assume here that this property remains true in the full theory with interaction. 
In |23| it has been argued that one could then expect to see a critical behaviour of 
the theory for r — > with r playing the role of the deviation of the temperature 
from the critical one, in the statistical physics of a system near a second order phase 
transition at T c , 

r~(T-T c )/T c . (5.29) 

We can then expect to see simple power behaviour of the correlation functions in X 4 
direction for r — > from general scaling and renormalization group arguments. In 



[23] various possibilities for the behaviour of the correlation length in X 4 direction 
are discussed. 

Let us assume here, as an example, that the r-dependent theory of section 4.3 
has a correlation length in X 4 direction similar to (|4.51|) 

U(r) ~ m- l (-r)-^ (5.30) 

for r — > 0, where m is a hadronic mass scale. We assume, furthermore, simple power 
behaviour of the correlation functions ( |4.7| ), ( |4.8| ) 

M- b (-iX 4 ,r) oc (-zX 4 M) a - 3 (-r)- 2 -H 

< a < 1, |e„| < 1, (5.31) 

for r — > and vrC 1 <C X 4 <ti \l^(r)\. The ansatz ( [5.31D corresponds to (see (A.l), 
(A.2)): 

M 2 (x°,r) ~ A l (x°M) a - 1 {-r)-^ £{ \ 

A' = const. (5.32) 



for r — > and m 1 .x° |^4( r )|- What are the consequences of ( 5.32|) for 



T 2 (ir], Q 2 ) of ( |3.36| ) and the structure function F 2 (xBj, Q 2 ) of ( |2.8| )? We have argued 
at the end of section 3.3 that the relevant region of the r integration in ( 3.36Q is for 
r w — f/2. Consider now the x° integral in ( |3.36| ) for this value of r. We have on 
the one hand the exponential damping factor exp(— x°/x(rj, r)), on the other hand, 
the matrix element A^^ , r) provides as cutoff the correlation length Z 4 (r). For 
r = —f/2 we get from fl3.34j ) and ( |5.30|) for 77 ^> Q: 

- ~ <£_ 

r ~ 2r/ 2 ' 
x(ry,-ir) w 

/ 4 (-|f) « -^r. (5.33) 
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This implies 



x{r],-\r) > h{-\r) 
x{r),-\f) < h{-\f) 



for Q < 2m, 
for Q > 2m. 



(5.34) 



Thus, for r] ^> Q ^> m, the x° integration in ( |3.36| ) is effectively cut off at a; = 
x(r],—f/2) and the x° integral in (|3.36| ) will get its leading contribution from the 
scaling regime of A^^ , r) (|5.32|) . Inserting ( |5.32| ) in (|3.36| ) we get for 77 ^> Q ^> m 



iMQ 2 



V 



dr 



c 



7-0 



dx°(x ) 2 

x exp{-x°/x{m r))A (x° M) a -\-r) 

O 2 / M\ 3+a r 
%-( — ) T(2 + a)A' / dr(-r)-* £0 (r + 

M 6 \ T] J Jc 

47r^ r(l + a+| £o ) f q 2 ^ 



-2-a 



-a- n£0 



M 



r ±e„) 



\2M 2 y 



M 



— l+a+eo 



(5.35) 



The analytic continuation for arbitrary large \v\ is obtained by the replacement in 



V 



(5.36) 



For real positive v this leads to 



vW 2 {v,Q 2 ) = F 2 (x Bj ,Q 2 ) 

~ 47rAT(l + a+i£o)(r(|£o))" 1 

l2M 2 J 



x sin 



1(1 



a - e 



-a— £0 



(5.37) 



This is an interesting result. The simple scaling assumption for the matrix el- 
ement M.2 (|5.32|) as suggested by the analogy to critical phenomena leads to the 
behaviour for the structure function F 2 found in the Regge fit (see (2) and (4a) of 
|p2| ) for large Q 2 and small x B f 



F 2 (x B] ,Q 2 )~X (Q 2 ) 



2\l+e 



{XBjY 



•eo 



Here X , Qq, e are constants with 



e « 0.44, 



(5.38) 



(5.39) 



where 1 + Eq is the intercept of the hard pomeron. Setting a = in (|5.37 ) we get 
indeed the powers of Q 2 and x B j of ( |5.38| ). Comparing (|5.28| ) and ( |5.31|) we see that 
\eq can be interpreted as the anomalous part of a critical index in the language of 
statistical physics. 
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6 Conclusions 



In this article we have developed an approach for the theoretical description of the 
structure functions at small XBj which should allow truely nonperturbative calcu- 
lations for all Q 2 > 0. Of course, one can also perform perturbative calculations in 
this framework. We have introduced an effective, r-dependent theory (see sect. 4) 
in Minkowski and Euclidean space, starting from QCD in Feynman-'t Hooft gauge. 
We have argued that the small xgj behaviour of the structure functions is related 
to the small r behaviour of the effective theory and that the limit r — > corre- 
sponds to a critical point. In the vicinity of this critical point we can expect to 
see power behaviour of the relevant matrix elements Aii^{x° ,r) (|3.28 ) with certain 
critical indices. We have shown that this leads to Regge behaviour of the structure 
functions at small XBj and large Q 2 , as observed in the phenomeno logical fits of 



[22]. In this way the intercept of the hard pomeron is related to a critical index 
of the r-dependent theory. We emphasize that in principle both the hard and the 
soft pomeron contributions to the structure functions should be calculable in our 
approach. 

The idea that the small XBj behaviour of the structure functions may be related 
to some kind of critical behaviour has been suggested by various authors. Our present 



article follows the ideas presented in [23]. The critical point of our effective theory 



for r — > 0, if it is indeed confirmed, would be completely analogous to a point of a 
second order phase transition in statistical physics. Quite a different type of critical 
behaviour, self-organised criticality, was suggested for the small xbj behaviour of the 
structure functions in | 31| . Criticality of the photon wave function in connection with 
a dipole model was suggested in It is also interesting to note that perturbative 
calculations in the leading logarithmic approximation, that is in the framework of 
the BFKL equation |H|, lead to conformal invariant structures for the amplitudes 
[32, 33 1 . One can suspect that this conformal invariance could have its origin in some 
sort of critical behaviour of an effective theory. 

Coming back to our present article we have presented here also our results on 
some more technical issues. 

We have discussed various ways of analytic continuation of the Compton am- 
plitude in the v and r planes. We found the representations ( p.29j )-( j3~36l ) the most 
useful ones. The r-dependent theory in Euclidean space was studied in particular 
in the Feynmann-'t Hooft gauge. The final goal of our approach is to make a truely 
nonperturbative calculation of the matrix elements (|4.26| ), ( |4.27|) in the r-dependent 
Euclidean theory. This could be based for instance on exact renormalization group 
methods (see |54, 35] for a review) or on lattice methods. Certainly, it will not be 



an easy task, since our Euclidean action ( |4.29| ) has an imaginary part for r ^ 1. 
But suppose that one can indeed study in this way the theory first for real r with 
< r < 2, deduce critical behaviour for small r and establish scaling relations as in 



(|5.31| ). Making then the analytic continuation in r and putting everything into the 
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theoretical machinery developed in this paper would mean a calculation of the small 
XBj and large Q 2 behaviour of the structure functions. Both the functional depen- 
dences and the absolute normalisation of the structure functions could be obtained 
in this way. Of course, much remains to be done. 
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Appendix A: Relation of M a i and -M 12 

From the definitions ( |3.16| ) and (|4.1|), ( |4.2| ) and using ( |3.5| ) we get after some straight- 
forward but lengthy algebra: 

U-( o n r o & 6 ( 1 " r ) ° 2 
M a (x°,r) = {-3- 



(9a; ) 2 x° dx°dr 
3r(2-r) d 2 6r(2 - r) 9 



(x ) 2 (<9r) 2 (x°) 2 (l - r) <9r 



}Mi(iV) 



r o.rf 52 6(1 -r)M 2 d 2 3r(2 - r) - 2 „ r2 <9 2 



(dx ) 2 x° dx°dr (x ) 2 (dr] 



2 



M^(x°,r) = -{^Y 2 ^ - J^^}{Ml(x ,r) + M 2 M^(x°,r)}. (A.2) 

Appendix B: Effective Lagrangian density in co- 
variant gauges 



Here we give the details of the derivation of C r (x) ( 4.14Q . From the original La- 



grangian £ (|4.9| ) we get the canonical momenta in the standard way. It is convenient 
to write first C separating time and space components of fields. 

C = l(G aj + «9 i G a0 -gf abc G b0 G cj ) 
(G aj + d j G° a -gf aV& G ,/0 G^) 

-^(G a0 + 9 j G aj )(G' a0 + d k G ak ) 

+ E{9|T°9-^7 9 
g 

+q( 1 2 i> dj -trfGfe - m q )q} 

+} a <j) a - (d^ a )d j( p a 

-gfa b c{4> a G b0 <p c + {djr)G b ^ c }. (B.l) 



We get now easily 



BC 1 • 

n ^ = ^ = -^ a0 + ^'), (B-2) 

U ° aJ = dG^ = ^ + d]Ga ° ~ 9fabcGb ° GCj > (R3) 
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fir 

n £a = ^ = ^-gf abc G b0 ct> c . 



Solving ( p.2|) , (|B.3|) , (|B.5| ) for the time derivatives of the fields, we 

G a0 = -£II Ga o - djG aj , 

G aj = U Ga3 - d 3 G a0 + gf ab cG b0 G c \ 

} a = IV, 

The Hamiltonian density is 

n = u Ga0 G a0 + n G aiG aj + 53(n 9 g + ^n ? ) 

or written out explicitly: 

IH. = —^UcaoIlQaO — YlQaodjG a:> 

+lU Ga jIL G aj + \G a jk G ajk 
+U Ga] (-d 1 G a0 + gf abc G b0 G^) 

+ E k-f 7 j dj +9YG a ^ + m q )q\ 

q > 

+n r n^ + ^ g f ahc G m <\> c 

+(d j r)d J r+gfaUd 3 r)G bj <p c . 

The — 3 component of the canonical energy-momentum tensor is 

T 03 = U Ga0 d 3 G a0 + U Ga ,d 3 G aj 
+U q d 3 q + (d 3 q)U q 

+n r 9 3 a + (9 3 a )n^. 

We have, furthermore, 



P 3 = / d 3 xT 03 = - / d 3 x% 



03 
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and from the definitions ( |4.4[ ), ( f4.12p 

H r = P° - (1 - r)P 3 



d 3 x{^ + (i-r)r 03 } 



(B.H) 
(B.12) 



Thus, we get 

H r = n+(l-r)T 03 . 

Inserting here ( P-8| ) and (|B.9|) leads to (|4.13|) . 

To prove that Ji r (4.13) is the Hamiltonian density to the Lagrangian C r (x) 
(|4.14j ) we proceed again in the standard way. We first write down C r (x) separating 
time and space components of fields. 

C r = l(G a J +d3G a0 -gf abc G b0 G c J -(l-r)d 3 G a J) 
(G aj + djG" - gfab>ciG b ' G c ' ] - (1 - r)d 3 G a] ) 

-±{G a0 + djG ai - (1 - r)d 3 G a0 ) 

(G a0 + d k G ak - (1 - r)d 3 G a0 ) 

+ E{?57 ?-^7°9 
i 

+g(|7 J dj -(1 - r)| 7 ° d 3 "27^4- - m 9 )g} 
+ (0 a - (1 - r)9 3 a )(0 a - gf abc G m <p c - (1 - r)«9 3 a ) 
-{djFMF - gfaUdj^G^r- 

From ( |B.13| ) we find the new canonical momenta as 

dC r 



rigao 



dG a0 
_1 
i 



HG 00 + 0jG aj ' - (1 - r)<9 3 G a0 ) 



n 



G a 3 



n 9 
it 



8C r 



-gfabcGV -(l-r)d 3 G a \ 

-(l-r)d 3 r, 
-gf abc G h0 <p c -{l-r)d 3 <p a 



(B.13) 



(B.14) 



(B.15) 



(B.16) 
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and Tl q , Hq staying as in ( |B.4| ). 

Solving ( |B.14| ), ( |B.15|) , ( |B.16|) for the time derivatives of the fields, we get: 



G a0 
G aj 



-fn G „o - d,jG ai + (1 - r)d 3 G a0 , 



U Ga3 - d,G a0 + gf ahc G m G C3 + (1 - r)d 3 G a i 



60, 



(B.17) 



LV + (l-r)<9 3 a , 
n, a +^ afec G 6 V c +(l-r)a 3 a . 



With this we find that 7i r (|4.13|) is precisely given as 

tC r = IlgaoG a0 + TiQajG^ 

+ ^(n 9 g + ^n 5 ) 

q 



(B.18) 



(B.19) 



Appendix C: Effective Lagrangian density in the 
temporal gauge 

Here we give the r-dependent Lagrangian and Hamiltonian densities in the temporal 
gauge. We start from the standard Lagrangian of QCD and add a total divergence 



term (see appendix E of f25|1 ) 



C 



\Gl p G aXp 



+ ^2 ( l(. i l X D\ - m q )q 
<? 

+d x \^{G aX d p G ap - G ap d p G aX ) - £ gf 7 A g 



Dy 



d x + igG a x -X a . 



With the gauge condition 



G a (x) = 0, (a=l, 



we get 



C 



±G aj G aj - i( 
-gf abc G^G bk d,G ck 



(CI) 



(C.2) 
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-\g 2 fabcfarsG^G Ck G^G Sk 

+ E {^7°? - + q d s -m, + <n J G aj y ) <?} . (C.3) 
The canonical momenta are 

dC ■ ■ 

U Ga] = —=- = G aj (C.4) 
.90 

and U q , Ilg as in (p.4|) . This leads to the following Hamiltonian density and 0-3 
component of the canonical energy-momentum tensor: 

= —TlQajUQaj 

+l(d k G aj )d k G a i 
-\{djG a i)d k G ak 



+gfa bc G a 'G bk d 3 G ck 
+\g 2 fa bc fastG h 'G ck G^G tk 

+ E {<? (-^7 J dj +m q - g^G»^j g} , (C.5) 

T 03 = UQajd 3 G^+J2{q^ dsq]. (C.6) 

With this we get for the r-dependent Hamiltonian density 
H r = H+(l-r)T 03 

= ni 0) + n ln \ (c.7) 

Hi 0) = ^UQa.UQa, + (1 - r)~U.Qa 3 d 3 G aj 

+\{d ] G ak )d ] G ak 
~(djG aj )d k G ak 



+ E {? (-^7 j dj +(1 - r) 1 -! ds g} , (C.f 



H Int = gf abc G aj G bk djG ck 



+\g 2 f ab cfastG b 'G ck G^G tk 
-E{^7^yg}. 



(C.9) 
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In the standard way we find from (U.8), ( |C.9| ) the r-dependent Lagrangian density 



C r = 4 0) +£ Int , (C.IO) 
4 0) = \{G aj - (1 - r)d 3 G aj )(G aj - (1 - r)d 3 G aj ) 

~(d k G aj )d k G aj 

+ E {? Ha -(1 - r)^ 7 ° g 3 -^) g} , (C.ll) 

We also find with £ from (|C.3| ). 

£ r = £ - (1 - r)G aj d 3 G aj 

+ ^(l-r) 2 (d 3 G a ^)d 3 G^ 



(l-r)E{^7°S 3 g}. (C.13) 



Appendix D: Euclidean Lagrangian density in co- 
variant gauges 

Let H r be as in ( |4.12j ), to be precise, we set x° = 0: 

H r = f d 3 xH r (x)\ x o =0 . (D.l) 



We define the Euclidean field operators with 

V(X 4 ) = exp(X A H r ) (D.2) 



as follows: 



G a Ej {X) = -V(X 4 )G aj (x)V{-X 4 ), 

G a m {X) = -iV{X 4 )G a0 {x)V{-X 4 ), 

q E (X) = V(X 4 )q(x)V(-X 4 ), 

q E (X) = V{X 4 )q{x)V{-X 4 ), 

</>%{X) = V(X 4 )0 ffl (x)\/(-X 4 ), 

4>%(X) = V{X 4 )4> a {x)V(-X 4 ), 

n G a.(x) = -v(x 4 )u Ga ,(x)v(-x 4 ), 

n G%4 (x) = iV(x 4 )u Ga0 (x)v(-x 4 ), 

n 0|; (x) = v{x A )u^{x)v{-x 4 ), 

IL*>(X) = V(X 4 )U, a (x)V(-X 4 ). (D.3) 
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Here G a3 (x) etc. on the r.h.s. of ( p.3| ) are the Minkowskian field operators at the 
point 

x = (X,0) 



and X = (X,X 4 ). 

Defining the Euclidean electromagnetic current operator as 

J E ,(X) = ^ {Q q q E {X) lE ,q E (X)} , 



(DA) 



(D.5) 



where Q q are the quark charges, we get with x as in ( D.4 ) 

J E4 (X) = V(X 4 )J°(x)V(-X,), 
J Ej (X) = -iV(X 4 )J\x)V(-X 4 ). 

The nucleon field operators in the Euclidean theory are obtained from (|4.16|) : 



(D.6) 



A Es (p,r) 
A Es (p,r) 



V{t)A s { p ,0)V(-t), 
V(r)Al(p,0)V(-r), 



± 



in 



(D.7) 



where we always take p (|2.18 ) for a nucleon at rest. We have 

lim A f Es (p,r) | O)exp(-Mr) = \N(p,s)), 

r — > — oo 

lim (0 | A Es {p, r) exp(Mr) = (N(p, s) | 



(D.8) 



In the standard way we get now for the matrix elements Ai a 6 (— 1X4, r) defined 
in for X 4 > 



M-(-iX A ,r) 



-^2(N(p,s) | (-1)J^(0,X 4 )J^(0,0) | N{p,s)) 

* s 

° — »- + 00 

(0 | A s (p,r / )(-l)J^(0,X 4 )J^(0,0)At(p,r i ) | 0). (D.9) 



MTi-iX^r) 



-]T(iV(p, s ) | J S4 (0,X 4 )J S4 (0,0) I iV(p,s)) 

* s 



(0 I A a (p,r / )Ji ! 4(0,X4)J^(0,0)At(p,r i ) | 0). (D.10) 
The Euclidean Hamiltonian density is obtained from ( |4.13|) with x from ( p.4|) 

H Er (X) = V(Xt)-Hr(x)V(-X4). (D.ll) 
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Using here ( p.3| ) leads to 



7i, Er (X) = — £ Ha a + —Hqo Tla a 



-I 



Ei 

+n G * i (id j G a EA + igU c G b m G c Ej + (1 - r)a 3 G^.) 

~ l ~4 ^Ejk ^Ejk 



2 

+n^n^ + n^(ip/ateG^ + (i - r)d^ a E ) 

-gfa b c(dA E )G b Ej <P E , (D.12) 

where all fields on the r.h.s. are at point X and we define the Euclidean field strength 
tensor as 

G Efll/ — 9fiG Eu — d v G a Elx — gf a bcG Ep G Ev . (D.13) 

It is now straightforward to verify that 7i Er is the Euclidean Hamiltonian density 
to the Lagrangian density C Er (|4.22|) . Indeed, starting from ( |4.22|) we define the 
Euclidean canonical momenta as follows, where G E4 = d^G E ± etc. 

.d£ Er 



9G EA 



n 



|{44 + ^ r (l-r)¥U 



. dC Er 

l : = 11 G 



E.1 



iG E j ~ idjG E4 — ig fabcG b E4 G c Ej — (1 — r)d 3 G a Ej , 



dC Er 






d<P% 








= 4% 


-(l-r)d 3 4>%, 


dC Er 


= % 








— ig fabcG E4 (j) C E 


dC Er 
dq E 


= n 9£ 


i 


dC Er 
dq E 


= n fe 


i 

= —^lE^lE- 



(l-r)d 3 r E 



(D.14) 
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The Hamiltonian density ( |D.12| ) is now obtained from (4.22) and ( p.14 ) as 



1 

-m^J% + + C Er , (D.15) 

where the ^-derivatives Ge^ etc. have to be considered as functions of the canonical 
momenta and the fields by inverting ( p.!4| ). 



Having derived the connection of the Euclidean Hamilton and Lagrange densities 
( |D.12| ) and ( [4.22 ) we can use the standard procedures of the path integral formalism 



to show that the matrix elements (|D.9| ) and ( |D.10| ) can be represented by the path 
integrals (Q) and (|Q7| ). 
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Figures 




> Rev 



Figure 2: The z/-plane with the position of the nucleoli poles (x) and the cuts. The 
funtions Tf(u,Q 2 ) (j = 1,2) are obtained along the short dashed line, Tj Ct (u,Q 2 ) 
along the long dashed line. 
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Figure 3: Cut structure of M(x°, 
plane for < e « 1. 



r, e) (a) and M.*(x°, r*, e) (b) in the complex r- 
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Figure 4: Cut structure of Aij(x°,r) and the curve C in the complex r-plane. For r 



sec 
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